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Abstract. We derive the Bethe ansatz equations describing the complete 
spectrum of the transition matrix of the partially asymmetric exclusion process 
with the most general open boundary conditions. By analysing these equations 
in detail for the cases of totally asymmetric and symmetric diffusion, we calculate 
the finite-size scaling of the spectral gap, which characterizes the approach to 
stationarity at large times. In the totally asymmetric case we observe boundary 
induced crossovers between massive, diffusive and KPZ scaling regimes. We 
further study higher excitations, and demonstrate the absence of oscillatory 
behaviour at large times on the "coexistence line" , which separates the massive low 
and high density phases. In the maximum current phase, oscillations are present 
on the KPZ scale t oc L~^/^. While independent of the boundary parameters, the 
spectral gap as well as the oscillation frequency in the maximum current phase 
have different values compared to the totally asymmetric exclusion process with 
periodic boundary conditions. We discuss a possible interpretation of our results 
in terms of an effective domain wall theory. 
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The partially asymmetric simple exclusion process (PASEP) [1,2] is a model 
describing the asymmetric diffusion of hard-core particles along a one-dimcnsional 
chain with L sites. Over the last decade it has become one of the most studied 
models of non-equilibrium statistical mechanics, see [3,4] for recent reviews. This is 
due to its close relationship to growth phenomena and the KPZ equation [5], its use 
as a microscopic model for driven diffusive systems [6] and shock formation [7], its 
applicability to molecular diffusion in zeolites [8], biopolymers [9 11], traffic flow [12] 
and other one-dimcnsional complex systems [13]. 

At large times the PASEP exhibits a relaxation towards a non-equilibrium 
stationary state. An interesting feature of the PASEP is the presence of boundary 
induced phase transitions [14]. In particular, in an open system with two boundaries 
at which particles are injected and extracted with given rates, the bulk behaviour in 
the stationary state is strongly dependent on the injection and extraction rates. Over 
the last decade many stationary state properties of the PASEP with open boundaries 
have been determined exactly [3,4, 15, 16, 19-22]. 

On the other hand, much less is known about its dynamics. This is in contrast 
to the PASEP on a ring for which exact results using Bethc's ansatz have been 
available for a long time [23-25]. For open boundaries there have been several 
studies of dynamical properties based mainly on numerical and phenomenological 
methods [26 30]. Very recently a real-space rcnormalization group approach was 
introduced, which allows for the determination of the dynamical exponents [31]. 

In this work, elaborating on [32] , we employ Bethe's ansatz to obtain exact results 
for the approach to stationarity at large times in the PASEP with open boundaries. 
Upon varying the boundary rates, we find crossovers in massive regions, with dynamic 
exponents z = 0, and between massive and scaling regions with diffusive {z = 2) and 
KPZ {z = 3/2) behaviour. 



The dynamical rules of the PASEP are as follows. At any given time t each site is 
either occupied by a particle or empty and the system evolves subject to the following 
rules. In the bulk of the system (i = 2,...,L— l)a particle attempts to hop one site 
to the right with rate p and one site to the left with rate q. The hop is prohibited if 
the neighbouring site is occupied. On the first and last sites these rules are modified. 
If site z = 1 is empty, a particle may enter the system with rate a. If on the other hand 
site 1 is occupied by a particle, the latter will leave the system with rate 7. Similarly, 
at i = L particles are injected and extracted with rates S and /3 respectively. 

With every site i we associate a Boolean variable Tj, indicating whether a particle 
is present (rj = 1) or not (r^ = 0). Let |0) and |1) denote the standard basis vectors in 
C^. A state of the system at time t is then characterized by the probability distribution 
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Figure 1. Transition rates for the partially asymmetric exclusion process. 
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where 

L 

\r) = \n,...,TL)=^\ri). (0.2) 

i=l 

The time evolution of \P{t)) is governed by the aforementioned rules and as a result 
is subject to the master equation 

^^\P{t))=M\P{t)), (0.3) 

where the PASEP transition matrix M consists of two-body interactions only and is 
given by 

M = Y^ I®^-^ ®M® I^L-k-i + ^ j^L-i ^ J0L-1 ^ ^Q_4^ 

k 

Here, / is the identity matrix on and M : (g) ^ (g) is given by 

'0 


,0 

The boundary contributions mi and mi describe injection (extraction) of particles 
with rates a and d (7 and /3) at sites 1 and L. In addition, as a tool to compute 
current fluctuations [33,34], we introduce a fugacity e^ conjugate to the current on 
the first site. The boundary contributions then are 

Strictly speaking, with the inclusion of A the matrix M is no longer a transition 
matrix of a stochastic process. In the following however, we will still refer to M as 
the transition matrix of the PASEP also for nonzero values of A. 

At A = 0, the matrix M has a unique stationary state corresponding to eigenvalue 
zero. For positive rates, all other eigenvalues of M have non-positive real parts. 
The large time behaviour of the PASEP is dominated by the eigenstates of M with 
the largest real parts of the corresponding eigenvalues. In the next sections we will 
determine the eigenvalues of M using Bethe's ansatz. The latter reduces the problem 
of determining the spectrum of M to solving a system of coupled polynomial equations 
of degree 3L — 1. Using these equations, the spectrum of M can be studied numerically 
for very large L, and, as we will show, analytic results can be obtained in the limit 
L — > 00. 



1. Relation to the spin-1/2 Heisenberg XXZ chain with open boundaries 

It is well known that the transition matrix M is related to the Hamiltonian H of the 
open spin-1/2 XXZ quantum spin chain through a similarity transformation 

M = -^U-^HxxzU^, (1.1) 
where f/^ and Hxxz are given by (see e.g. [19,20]), 

^'^ =<S){1 (1-2) 

i=l ^ ^ 
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L-l 



xxz 



(1.3) 



+ Bi+Bl 



The parameters A and h, and the boundary terms Bi and Bl are related to the 
PASEP transition rates by 



Bl 



1 



Br = 



1 



(a + 7 + (a - 7)<Ti - 2ane^a^ - 2'yiJ,~^e~^a^) , 



(1.4) 



Here cr^, uj and ctj are the usual Pauli matrices, — (ct| ± ia^)/2, and /z is a free 
gauge parameter on which the spectrum does not depend. 

The expression in (1.3) is the Hamiltonian of the ferromagnetic Uq{SU{2)) 
invariant quantum spin chain [35] with added boundary terms Bi and Bl- We 
note that the boundary terms contain non-diagonal contributions {af, a^) with L- 
dependent coefficients. In the absence of the boundary terms the spectrum of the 
Hamiltonian is massive, i.e. there is a finite gap between the absolute ground states 
and the lowest excited states. As is shown below, the boundary terms lead to the 
emergence of several different phases. Some of these phases exhibit a spectral gap in 
the limit L — > oo in the sense that there is a finite gap in the real part of the eigenvalue 
of the transition matrix. Other phases feature gapless excited states. 

In order to make contact with the literature, we start by relating the PASEP 
parameters to those used in previous analyses of the spin- 1/2 XXZ quantum spin 
chain with open boundaries. For the latter we employ notations similar to [36], in 
which the XXZ Hamiltonian reads (up to a constant shift in energy) 



H 



Bl = 



1 

2 ^ 

i=l 



+ ^^^+1 - COST? (af <7^+i - 1)] + Bl + Bl, 



smr] 



cos w_ + cos ^_ 



Br 



siu ;/ 



COS a;+ + cos (5+ 



-(cosa;_ — COS (5_) erf 

+ cos6ia^ + sin^icr^ — sina;_ 
- ^ (cos w+ — COS (5+ ) (t£ 

+ COS 62 cr£ + sin 62 (j\ — sin a;+ 



(1.5) 



Comparing this with the Hamiltonian (1.3) arising from the PASEP, we are led to the 
following identifications 



-le" 



= -ie"^+ 



Exact Spectral Gaps of the Asymmetric Exclusion Process with Open Boundaries 5 



Purthermore, S± are determined via the equations 

faj sin 7^ I IJd smr) 



pq cos u)- + cos 6- y pq cos u)+ + cos S+ 



(1.7) 



As the spectrum of H cannot depend on the gauge parameter /x, it follows from (1.6) 
that the spectrum depends on 9i and 62 only via 9i — 92- 

Although it has been known for a long time that H is integrable [37, 38] , the 
off-diagonal boundary terms have so far precluded a determination of the spectrum of 
H by means of e.g. the algebraic Bethe ansatz [39]. However, recently a breakthrough 
was achieved [40-42] in the case where the parameters satisfy a certain constraint. In 
the above notations this constraint takes the form 

cos(6ii - 62) = cos((2fc + l)r) + uj-+ uj+). (1.8) 

Here k is an integer in the interval \k\ < L/2. In terms of the PASEP parameters the 
constraint reads 

(0^+2^ - e^)(a/je^ - Q^-^'^-^jS) = 0. (1.9) 

For given k and A the constraint (1.9) can be satisfied in two distinct ways. Either 
one can choose Q as a root of the equation (5^+^*^ = e^, or one can impose a relation 
on the boundary and bulk parameters such that the second factor in (1.9) vanishes. 
Curiously precisely under the latter conditions, the DEHP algebra [15] for the PASEP 
fails to produce the stationary state [20,43]. 

Either choice results in a constraint on the allowed rates in the PASEP. In 
particular this implies that it is not possible to calculate current fiuctuations for general 
values of the PASEP transition rates from the Bethe ansatz equations presented below. 

However, in the ease A = it is possible to satisfy the constraint (1.9) for arbitrary 
values of the PASEP parameters by choosing k = —L/2. As the condition (1.8) at 
first sight appears not to have any special significance in terms of the XXZ chain, it 
is rather remarkable that it can be identically satisfied for the PASEP. 



1.1. Symmetries 

The spectra of M and H are invariant under the particle- hole and left-right symmetries 
of the PASEP: 

• Particle-hole symmetry 

a 7, 13-^5, P q, A — > — A. (1-10) 

• Left-right symmetry 

a -i-^ 5, /3 7, P q, A ^ — A. (l-H) 

Both particle-hole and left-right symmetries leave the two factors in (1.9) invariant 
individually. There is a third symmetry which leaves (1.9) invariant, but interchanges 
the two factors 

• Gallavotti-Cohen symmetry [44,45] for the PASEP 
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The combination of (1.12) with a redefinition of the gauge parameter, 

fi^ l^e^Q-^+\l^, (1.13) 




corresponds to the interchange 6i 62- It is shown in the next section that if the 
constraint (1.8) is satisfied, the spectrum of H no longer depends on the difference 
9i — 02- Hence (1.12) is not only a symmetry of the constraint equation (1.8), but also 
of the spectrum of H. 

The Gallavotti-Cohen symmetry (1.12) implies the fluctuation theorem [46,47] 
for the probability Pl{J, t) to observe a current J on the first site at time t, 

We emphasize the L dependence in (1.12) and (1.14), which disappears for symmetric 
hopping, Q = 1 [48]. It is further clear from (1.12) that the current vanishes J = 
when the detailed balance condition 

is satisfied [43]. Precisely the same Gallavotti-Cohen symmetry as in (1.12) was 
observed for the zero-range process (ZRP) with open boundaries [49]. This model 
is equivalent to the PASEP on an infinite lattice but with a fixed number of i + 1 
particles and particle dependent hopping rates. 

As a final remark, we note that it was realized in [36] that the constraint 
(1.8) has in fact an algebraic meaning and corresponds to the non-semisimplicity of 
an underlying Temperlcy-Licb algebra, and implies the existence of indecomposable 
representations. Condition (1.8) can be interpreted as a generalized root of unity 
condition for this algebra, and it implies certain additional symmetries for the XXZ 
chain. It would be of interest to understand the impact of these symmetries on the 
spectrum of the PASEP with open boundaries^. 

2. Bethe ansatz for the XXZ Hamiltonian 

The first Bethe ansatz results pertaining to the spectrum of H (1.5) were reported 
in [40,41]. Subsequently it was noted on the basis of numerical computations in [42] 
and an analytical analysis for a special case in [50], that these initial results seemed 
incomplete. Instead of one set of Bethe ansatz equations, one generally needs two. 
Further developments were reported in [51,52]. 

In order to simplify the following discussion we introduce the notations 
2sin« sina;+ 

a± = ^. 2.1 

cos U}± + COS 0± 

When the constraint (1.8) is satisfied for some integer k, the eigenvalues of H can be 
divided into two groups, Ei{k) and E2{k) 

L/2-l-k . 2 



X However, we have checked that for the TASEP with open boundaries the spectrum consists of 
singlets only, confirming [64], and hence is diagonalisable. 
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L/2+k 



2 sin T] 



(2-3) 

^-^ cos 2vi — cos r? 

j=i J 

Here the complex numbers {ui} and {vi} are solutions of the coupled algebraic 
equations (i = 1, . . . , L — 1) 



where 



w{Vi) 



w{u) 



2L 



K_{-Ui - uj-)K+{-Ui - w+) J-j*- S{-Ui,'i 
K-{vi)K+{vi) ^jj'' S{vi,Vj) 



K_{-Vi)K+{-Vi) J-j^ S{-Vi,vj) 



(2.5) 



sin(ry/2 + Ui) 



S{u,v) = cos2v — cos(277 + 2u), 



sin(r//2 - «,) ' 

K±{u) = cos^i + cos(ry + uj± + 2u). (2.6) 

We will now rewrite these equations in terms of the PASEP parameters for which 
it is convenient to employ the notations 



z = — e 



2iu 



c 



^2 it) 



and 



a, 7 
.± 



p - g - a + 7, 



We then find that, 

• equations (2.2) and (2.4) are rewritten as 



^Ei{k) = a + P + -f + 5+ Yl 



- {Q-zj){Qzj-iy 



(2.7) 

(2.8) 
(2.9) 



(2.10) 



ZjQ - 1 



Ki{z,a,^) 



(2.11) 
(2.12) 



Q - 2,- 



Here i4ri(2;) = i^i(2, a, 7)ifi(z, /3, i5) and 

{z + Qk+^^){z + Qk-,^) 



{Ql<a,iZ + l)((5Ka,72 + 1) 



• the second set of equations, (2.3) and (2.5), becomes 
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Q-Cj 



2L 



L/2+fe 



^2(0)= n 



where i^2(C) = ^2(C,a,7)^2(C,/3,^) and 

, ('vT,-,( + Q)('vT.-A+Q) 

K2{z,a,j) 



.L-1, 
(2.14) 

(2.15) 



(QC + Ki,7)('9C + '«a,7)' 

As wc have mentioned before, in the case of the PASEP, and for A = 0, the 
constraint (1.9) can be satisfied by either considering symmetric hopping, Q = 1, or 
by choosing k = —L/2. 

3. Bethe ansatz equations for the "generic" PASEP 

Inspection of the second set of equations (2.13) for the choice k = —L/2 reveals that 
there exists an isolated level with energy Sq = 0. This is readily identified with the 
stationary state of the PASEP. Furthermore, all other eigenvalues £ of the transition 
matrix M follow from the first set of equations (2.10) and (2.11), and are given by 

L-1 /^T ,\2 



£ = -a-(3-^-d-Y^ 



p{Q'-l) 



^ {Q-zj)iQzj-iy 
where the complex numbers zj satisfy the Bethe ansatz equations 



(3.1) 



ZjQ - 1 



Q - Zj 



2L 



L-1 



Kiz,) = n ?r 



Zl ZjZli 



- 1 



Zj - ziQ^ ZjZi - Q"^ 



1...L-1. 



Here K{z) = K{z, a, j)K{z, p, 6) and 



Ki{z,a,'y) 



{z + Qk+){z + ^ 



(3.2) 



(3.3) 



{QkZi^^Z + l){QKa,-fZ + 1) 

In order to ease notations we set from now on, without loss of generality, p = 1 and 
hence Q = sjq. 

4. Symmetric Exclusion Process (SEP) 

The limit of symmetric diffusion is quite special and we turn to it next. We can 
obtain this limit either by taking Q ^ \ and leaving k unspecified in the equations 
(2.10), (2.11) and (2.13), (2.14), or by setting k = -L/2 and studying the limit of the 
equations (3.1), (3.2) for the generic PASEP. The two procedures lead to the same 
results and we will follow the second for the time being. 

Taking the limit Q — > 1 in the equations (3.1), (3.2) for the general PASEP we 
observe that 

ZjQ^ 



zQ-1 
Q-z 
ZjZiQ"^ ■ 



-1 



Zl 



Zi 



1 



ZjZl - 



Q^zi 

K{z,a,'j) 



az + ^ 
'yz + a' 



(4.1) 
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We conclude that in this hmit the spectral parameters z must fulfill 

az^^z^ (4.2) 

jz + a oz + p 

It is easy to see that the only solutions to this equation are z = ±1. Having established 
the leading behaviour of the roots of the Bethe ansatz equations in the limit Q ^ 1 
we now parametrize 

Zj=±l + iXj{Q^-l) , (4.3) 

then substitute (4.3) back into the Bcthc ansatz equations (3.2), (3.1) and finally take 
the limit Q ^ 1. Choosing the plus sign in (4.3), we obtain with c{x) = {x + 2)/2x 

2L 



A, - i/2 
Xj + i/2 



Xj + i c{a + 7) 



L-l 

= n 



Xj — ic{a + 7) 
Xj - A; - 



Xj +ic{P + 6) 



A, 



Xj - A; + i 



Xj + Xi 



A,- + A; + i 



ic{P + 6)_ 



£ = —a — /? — 7 — i5 - 



1/4- 



(4.4) 



(4.5) 



On the other hand, choosing the minus sign in (4.3) we arrive at 



£ = —a — /3 — 7 



Ai + i Aj + A; + i 
5 . 



(4.6) 



(4.7) 



We observe that for this choice of sign in (4.3) we obtain only a single energy level. 

For both choices there is a subtlety: we have implicitly assumed that the Aj's 
remain finite when we take the limit Q ^ 1 and this need not be the case. A numerical 
analysis of the Bethe ansatz equations (3.2) for Q « 1 shows that we have to allow 
one or several Aj's to be strictly infinite, in which case they are taken to drop out of 
(4.4) §. This then leads to the following set of Bethe ansatz equations, in which we 
only allow solutions where all spectral parameters Xj are finite 

2L 



A, 



i/2V Xj+ic{a + j)Xj+ic{p + 6) 



Xj+i/2j Xj -icia + 'y) Xj -ic{l3 + S) 



N 



-j-j- Xj - Xi-iXj + A/ 



Xj - A; + i Aj + A; + i 



,iV. 



N 



-a- P - 5 



^A,^ + l/4- 



(4.E 



(4.9) 



Here A'' is allowed to take the values 1, 2, ... L — 1. 

Curiously, the limit of symmetric exclusion can be described by a second set 
of Bcthc ansatz equations. Taking the limit Q ^ 1 of (2.13), (2.14) and leaving k 



§ We have verified this prescription by solving the Bethe ansatz equations numerically for small 
systems in the vicinity of the limit Q — > 1. 
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unspecified we arrive at £ = or 

2L 



- i/2 \ Aj + i d{a + ) Xj + i d{ J + 5) 



\j+i/2J Xj-id{a + 'y)Xj-id{l3 + 6) 



N 

n 

1^3 



Aj - A; - i Aj + A; - i 
Aj - Ai + i A, ■ + A; + i 



N 



^A2 + l/4- 



Here 



d(a;) — c(— x) 



2a; 



(4.10) 



(4.11) 



(4.12) 



Apart from a constant shift in energy the two sets of equations are related by the 
simultaneous interchange a —7 and (3 <-> —5. Importantly, equations (4.10) 
coincide with the Bethe ansatz equations derived in [53] by completely different means. 
Numerical studies of small systems suggest that either set (4.4) or (4.8) gives the 
complete spectrum of the Hamiltonian. 

Interestingly, the Bethe equations (4.8) and (minus) the expression for the energy 
(4.9) are identical to the ones for the open Heisenberg chain with boundary magnetic 
fields [54,55] || 



H : 



j+1 - 7 



1 



1 

i7 



(4.13) 



-2E 

where ^_ = — 1/(0 + 7) C+ = + if the reference state in the Bethe ansatz 

is chosen to be the ferromagnetic state with all spins up. On the other hand, the Bethe 
equations (4.10) and (minus) the expression for the energy (4.11) are obtained when 
the reference state is chosen as the ferromagnetic state with all spins down. This on 
the one hand shows the equivalence of the two sets of Bethe equations and on the other 
hand establishes the fact, that the spectrum of the SEP with open, particle number 
non-conserving boundary conditions (which corresponds to the z-component of total 
spin not being a good quantum mimbca' in the spin-chain language) is identical to 
the spectrum of the open Heisenberg chain with boundary magnetic fields, for which 
5^ is a conserved quantity. This spectral equivalence is more easily established by 
means of a similarity transformation, see section 6.5.1 [4] and [56] for the case at 
hand and [57] for a general discussion on spectral equivalences for conserving and 
non- conserving spin chains. 

The first excited state for the SEP occurs in the sector A'' = 1 of (4.10). The 
solution to the Bethe ansatz equations for large L is 

TT 

I2Z 



Ai 



- - - [d{a + 7) + d{[3 + 6)] 

TT TT 



+ 7) - 2(i3(a + 7) + d{(3 + 5)- 2d^(f3 + S)] + . .(.4.14) 



The corresponding eigenvalue of the transition matrix scales like L ^ with a coefficient 
that is independent of the boundary rates 

S,{L) = -^+OiL-% (4.15) 



II We note that for symmetric diffusion the eigenvalues of the transition matrix are simply minus the 
corresponding eigenvalues of the Heisenberg Hamiltonian. 
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We conclude that for the SEP the large time relaxational behaviour is diffusive and 
universal. 

5. Totally Asymmetric Exclusion Process (TASEP) 

We now turn to the limit of totally asymmetric exclusion Q = and set 7 = (5 = in 
order to simplify the analysis. 

5.1. Stationary State Phase Diagram 

The stationary state phase diagram for the TASEP was determined in [15,16,58] and 
features four distinct phases, see Figure 2: 

1. A low density phase for a < 1/2, a < /3. 

In the low-density phase the current in the thermodynamic limit is equal to 
J = a{l — a) and the density profile in the bulk is constant p = a. 

2. A high density phase for /3 < 1/2, /3 < a. 

Here the current in the thermodynamic limit is equal to J = f3{l — [3) and the 
density profile in the bulk is constant p ~ \ — (3. 

3. A coexistence line at (3 — a < 1/2. 

On the coexistence line the current is equal to J = a(l — a), but the density 
profile increases linearly in the bulk p{x) = a + (1 — 2a)x. 

4. A maximal current phase at a,/? > 1/2. 

This phase is characterized by the current taking the maximal possible value 
J = 1/4 and the bulk density being constant and equal to p — 1/2. 

In [16] a further subdivision of the low and high density phases was proposed on 
the basis of differences in the behaviour of the density profile in the vicinity of the 
boundaries. In the high-density phase this distinction corresponds to the parameter 
regimes a < 1/2 and a > 1/2, respectively. 




Figure 2. Stationary pliasc diagram determined by tlie current of the TASEP. 
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5.2. Analysis of the Bethe ansatz equations 

In order to determine the exact value of the spectral gap we will now analyse (3.1) 
and (3.2) in the limit L — > oo. 

After rescaling z Qz and setting j = 6 = 0, the Q ^ limit of equations (3.1) 
and (3.2) reads 

£ = -a-0-J2-^, (5.1) 
1=1 ^' ^ 



(5.2) 



In order to ease notations in what follows, we introduce 



5W =1M (7^). (5-3) 
gh{z) = In ( ) + In (^; + a) + In {z + b), (5.4) 



where 



a=--l, b=^-l. (5.5) 

a p 

The central object of our analysis is the "counting function" [59-61], 

1 1 ^"^ 

iYl{z) = g{z) + -gi,{z) + - ^ K{zi, z), (5.6) 

1=1 

where K{w,z) is given by 

K{w, z) = -~\nw + ln{l — wz). (5.7) 

Using the counting function, the Bethe ansatz equations (5.2) can be cast in 
logarithmic form as 

yL{zj) = ^Ij, j = l,...,L-l. (5.8) 

Here Ij are integer numbers. The eigenvalues (5.1) of the transition matrix can be 
expressed in terms of the counting function as 

£ = -a-(3-L\im (^iYHz) - g'iz) - igU^)) • (5-9) 

Each set of integers {Ij\j = I, . . . ,L — 1} in (5.8) specifies a particular excited 
state. In order to determine which set corresponds to the first excited state, we have 
calculated the eigenvalues of the transition matrix numerically for small systems of 
up to L = 14 sites for many different values of a and /?. By comparing these with 
the results of a numerical solution of the Bethe ansatz equations, we arrive at the 
conclusion that the first excited state always corresponds to the same set of integers 

Ij = -L/2 + j for j = l,...,L-l. (5.10) 

The corresponding roots lie on a simple curve in the complex plane, which approaches 
a closed contour as L ^ oo. The latter fact is more easily appreciated by considering 
the locus of reciprocal roots z^^ rather than the locus of roots Zj. In Figure 3 we 
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present results for a = /3 = 0.3 and a = /? = 0.7 respectively. The limiting shape of 

the curve is that of the cardioid, which can be seen as follows. Assuming that the last 
term in (5.6) is approximately constant as L ^ oo and using (5.8), (5.10) we find that 

exp(-g(z,)) = (zj/^ - zj'/' f = xe"''''^/^ {L ^ ^). (5.11) 

Here x = x{(^^P) is some constant. Parametrising z = pe'^ and multiplying (5.11) by 
its complex conjugate, we conclude that the roots lie on the curve defined by 

p2-(x + 2cos^)p+l = 0, (5.12) 

the defining equation of the cardioid. 




The shape of the locus of inverse roots depends on the rates a and /?. For example, 
in the case a = /3 = 0.7 shown in Figure 3(b), a cusp is seen to develop at the intercept 
of the curve with the negative real axis (which occurs at 2; = —1 in the limit L — > 00). 
In terms of the parameter x characterizing the cardioid the cusp develops at x = 4. 
In contrast, no cusp occurs for a = /? = 0.3 as shown in Figure 3(a). ^ We will see 
that this difference in the shape of the loci of inverse roots is reflected in a profoundly 
different finite-size scaling behaviour of the corresponding spectral gaps. 

In order to compute the exact large L asymptotics of the spectral gap, we derive 
an integro- differential equation for the counting function Yl{z) in the limit i — > 00. 
As a simple consequence of the residue theorem we can write 

where C = Ci + C2 is a contour enclosing all the roots Zj, Ci being the "interior" 
and C2 the "exterior" part, see Figure 4. The contours Ci and C2 intersect in 
appropriately chosen points ^ and It is convenient to fix the end points ^ and 

^* by the requirement 

Using the fact that integration from ^* to ^ over the contour formed by the roots is 
equal to half that over C2 — Ci we find, 

iYUz) = g{z) + y^{z) + K{w, z)Yl{w)dw 

1 f K{w,z)Yl{w) 1 f Kiw,z)Yl{w) 



*l This corresponds to having x > 4. 
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where we have chosen the branch cut of K{w, z) to he along the negative real axis. 

Im z 




Rez 



Figure 4. Sketch of the contour of integration C in (5.13). Tlie open dots 
correspond to the roots Zj and £, is chosen close to 2i_i and avoiding poles of 
cot{LYL{z)/2). 

Our strategy is to solve the integro- differential equation (5.15) by iteration and 
then use the result to obtain the eigenvalue of the transition matrix from equation 
(5.9). 

6. Low and High Density Phaises 

In the low and high density phases the locations of the end points ^ and ^* are such that 
a straightforward expansion of (5.15) in inverse powers of L is possible (see e.g. [62,63] 
and Appendix A). The result is 

\Yl{z) = g{z) + ^g^{z) + ^ j^^^ K{w, z)Yi{w)dw 



12L2 V y/(e) n(0 



= 9iz) + ^9h{z) + ^ f° K{w, z)Yi{w)dw 

+ i- K{w,z)Y[{w)dw + ^ J[ K{w,z)Yi{w)d 



12L2 V n(r) Y[{o J ' " ^^-^^ 

where the derivatives of K are with respect to the first argument. Wc note that here 
we have imphcitly assumed that y^(^) is nonzero and of order ©(L^). In order to find 
the eigenvalue of the transition matrix (5.9) up to second order in inverse powers of L, 
we will need to solve (6.1) perturbatively to third order. Substituting the expansions 



(6.2) 



n=0 



back into (6.1) yields a hierarchy of equations for the functions yn{z) of the type 



1 

Vniz) = gn{z) + ^ / _ ^(^' ^)yn{w) dw, 



(6.3) 
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where = Zc±iO. The integral is along the closed contour following the locus of the 
roots. The first few driving terms gn{z) are given by 

go{z) = -ig{z), 

9i{z) = -igh{z) + Ki + XiK(zc,z), ^g^^ 

52(2) = K2 + X2K{Zc, Z) + /i2if'(Zc, z), 

gsiz) = K3 + hK{zc, z) + fi3K'{zc, z) + u^K"{zc, z). 

We recall that the functions g and gb are defined in (5.3) and (5.4) and K{zc, z) = 
\n{z — z~^). The coefficients k„, A„, ^„ and are given in terms of ry„ defined by 
(6.2) as well as derivatives of y„ evaluated at z^- Explicit expressions are presented in 
Appendix B. 

These coefficients as well determined self-consistently by solving (6.3) 

and then imposing the boundary conditions (5.14). 



6.1. Small values of a and (3 



When a and /3 are small we assume that the singularities in gh{z), i.e. the points 
—a — 1 — l/a and —6 = 1 — 1/(3, lie outside the contour of integration. From the 
distribution of the reciprocal roots, Figure 3(a), we further infer that for small values 
of a and /3 the roots lie inside the unit circle. In particular we assume that Zc ^ —1 
and that the points dzl lie outside the contour of integration. 

We now proceed to solve (6.3) and then verify a posteriori that the above 
assumptions hold. 

The equation (6.3) for n = is solved by the simple ansatz 



yQ{z) = Ko+goiz) . 



(6.5) 



Substituting the ansatz into the integro-differential equation (6.3) for n = we find 
that 



Ko 



{—\nw + ln(l — wz)) 



w 



1 



dw 



= - i (- In(-^c) + 2 ln(l - 0e)) • (6.6) 

This in turn implies that the zeroth order term in the expansion of the counting 
function is given by 

2' 



yo{z) 



iln 



z 

Zc 



1 



(6.7) 



In order to derive this result we have made use of the following simple but useful 
identity (C denotes the contour of integration from z~ to z^) 

Inti; , f ln(l + 2;c/a;) if - a; outside C, 
ln(— a; — Zc) if — x inside C. 



27ri J,- 



dw 



(6.8) 



The integro-differential equations for n 
with the results 



1,2,3 are solved in an analogous manner, 



yi{z) = - iln 



z 1 



Zcl- z'^ 



Z- Zc 



-iAi 



z + a z + b 
Zc + azc + b 



+ Ki-nn{ab{-zc)-'^') 



(6.9) 
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»W = A.l„(^V§f-^-— ^ 



Z Zq / Z^ y z Zq Zq Zq 

+ K2-A2ln(-0e)-— , (6.10) 

Zc 

2/3W = A3ln(^) + (;^-^) 

Vz-ZcV \Zc zlJ\z-Zc^ Zc-Zc^J 

Ks - X3ln{-Zc) 



V3 / Z Zc 

z'c\{z-zc'r {zc-zc'r 

Zq Z^ 



(6.11) 



The coefficients k„, A„, fin and Vn are given in terms of Sn, rjn and derivatives of 
Un at Zc, see Appendix B and we now proceed to determine them. Substituting the 
expansions (6.2) into equation (5.14), which fixes the endpoints ^ and we obtain a 
hierarchy of conditions for yn{zc), e.g. 

Yl{0 = yoiO + jyiiO + ;^2/2(0 + • • • 

= yo{zc) + 2 iViizc) + y'o{zc){Si + irji)] + ... 

= ^-j- (6.12) 
Solving these order by order we obtain 
Ai = 2i, 

A3 = M2 = A2 = Kl = 0, 

1/3 = ZcHz = Zc{Zc - 1)^K2 = -iTT^-^c I , i 

with Lmdctcrmincd. Furthermore, the intersect of the solution curve with the 
negative real axis in the limit L ^ 00 is found to be 

1 



1X2 



(6.13) 



- ^. (6.14) 

Having determined the counting function, we arc now in a position to evaluate the 
corresponding eigenvalue of the transition matrix from equation (5.9) 

(X) = - a - - ^ - ^ + O (X-) 

i- Zc Zc — Zc ^ 

= ~a-(i+ ^ - -—^^—1- + O (L-3) . (6.15) 

l + ^^6 Vab-l/VabL^ \ j y j 

We conclude that to leading order in L, the eigenvalue of the transition matrix with 
the second largest real part is a nonzero constant. This implies an exponentially fast 
relaxation to the stationary state at large times. We note that due to the symmetry of 
the root distribution corresponding to (5.10) under complex conjugation £1 is in fact 
real. The domain of validity of (6.15) is determined by the initial assumption that —a 
and —b lie outside the contour of integration, i.e. 

—a < Zc and — b < Zc. (6.16) 



Exact Spectral Gaps of the Asymmetric Exclusion Process with Open Boundaries 17 



The parameter regime in a and /3 in which (6.15) is valid is therefore bounded by the 
curves 

-1 



/3c 



1 



a 



1-a 



1-/3 



1/3- 



1/3 



, 0<a<^, 



o<p<-. 



(6.17) 



6.1.1. Coexistence Line On the Hne f3 = a, the leading term in (6.15) vanishes and 
the eigenvalue with the largest real part different from zero is therefore 

7r^a(l - a) 1 



I -2a L2 



(6.18) 



This equation holds for fixed a < 1/2 and L oo. We note that there is a divergence 
for a ^ 1/2, signaling the presence of a phase transition. The inverse proportionality 
of the eigenvalue (6.18) to the square of the system size implies a dynamic exponent 
z = 2, which suggests that the dominant relaxation at large times is governed by 
diffusive behaviour. As shown in [29, 30] the diffusive nature of the relaxational 
mechanism can in fact be understood in terms of the unbiased random walk behaviour 
of a shock (domain wall between a low and high density region). Our results (6.15) 
and (6.18) for the phase domain given by (6.17) and the coexistence line agree with 
the relaxation time calculated in the framework of a domain wall theory (DWT) model 
in [30]. As we will show next, this is in contrast to the massive phases beyond the 
phase boundaries (6.17), where the exact result will differ from the DWT prediction. 
In Section 9 we present a modified DWT framework that allows us to understand the 
phase boundaries (6.17). 



6.2. Massive phase II: 

In this section we will treat the case — fe > Zc with —a < Zc as before. The case 
—a > Zc with —b < Zc is obtained by the interchange a ^ /3 in the relevant formulas. 

We need to solve the same integral equation (6.3) as before; in particular the 
driving terms g„ defined in (6.4) remain unchanged. However, when iterating the 
driving term gi{z) there is an extra contribution because the pole at —b now lies 
inside the integration contour, see (6.8). As a result the solution for yi{z) is now of 
the form 

'"^^ z + a z + b z + 1/b 
Zc + a Zc + b Zc + 1/b 

+ Ki-iln(a(-Zc)^'^i). (6.19) 

The forms of the solutions for for n > 2 remain unchanged. The coefficients of 
the various terms are again fixed by imposing the boundary condition (5.14), with the 
result 

Ai = 3i, 

As = A«2 = A2 = Ki = 0, 
3 

1^3 = ZcfJ-3 = o^c(^c - 1)^K2 



yi{z) 



iln 



z 1 



Zc 1 



, (6.20) 
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and ^3 again remains undetermined. We note that the difference in the functional 

form of 2/1(2;) compared to (6.9) affects the coefficients in all subleading contributions 
2/2(2)) 2/3(2;) etc. The intersect of the locus of roots with the negative real axis in the 
limit i ^ 00 is now given by 

-,1/3 



/abc 



-1/3 



a 



1 



a 



where we have defined 



bc = 



I -(3c 



(6.21) 



(6.22) 



The eigenvalue of the transition matrix with the largest nonzero real part is again 
determined from (5.9) 



l+2Ze 



1- Zc 

= - a - /3c + ■ 



L2 

47r2 



1 



1 + ai/3 aVs _ a-i/3 2,2 
The result (6.23) is valid in the regime 

0<a<l/2, /3c </?<!, 



+ O . (6.23) 



(6.24) 

and is seen to be independent of (3. Hence in this phase the relaxation to the stationary 
state at large times is independent of the extraction rate at the right-hand boundary 
of the system. 



7. MELximum Current Phase 



In the maximum current phase a, /? > 1/2 the above analysis of the intcgro-diffcrcntial 
equation (5.15) breaks down. The primary reason for this is that the locus of roots 
now closes at 2;c = — 1 for L — > 00 and this precludes a Taylor expansion of the kernel 
K{w,z) around w = ^. A more detailed discussion of the complications arising from 
2;c = — 1 is presented in Appendix A. 

In order to determine the large-L behaviour of the eigenvalues of the transition 
matrix with the largest real parts in the maximum current phase, we have resorted 
to a direct numerical solution of the Bethe ansatz equations (3.2) for lattices of up 
to L = 2600 sites. In order to facilitate a finite-size scaling analysis it is necessary to 
work with quadruple precision (32 digits in C) at large L. Wc find that the leading 
behaviour of the spectral gap in the limit L ^ 00 in the maximum current phase is 
independent of the boundary rates a and /3. 



7.1. Leading behaviour 

In Figure 5 we plot the numerical results for eigenvalue £1 of the first excited state 
for a = /3 = 0.7 as a function of the inverse system size on a double-logarithmic 
scale. The almost straight line suggests an algebraic behaviour at large L 

£i{L) ~ -eL-' {L 00). (7.1) 

A simple least-square fit of the graph in Figure 5 in the region 2200 < L < 2600 to a 
straight line gives a slope of s « 1.493, which is close to the expected value 3/2. 
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Figure 5. Double logarithmic plot of —£i{L) as a function of 1/L for a = j3 = 0.7 

A better result is obtained by extrapolating our finite-size data as follows. We 
first divide the data set for £i{L) into bins containing 11 data points each. The fc'^ 
bin Bk is defined by taking 20k < L < 20{k + 1). Within each bin we fit the data for 
£i{L) to a functional form 

£i{L)^-ekL-"'. (7.2) 

We have implemented this procedure for 100 < k < 128 and obtained a sequence of 
exponents Sk- We observe that the following least-square fit of the sequence Sk to a 
polynomial gives excellent agreement 

Sk w 1.4999949 - 0.7822533fc~i + 3.8014387fc~2. (7.3) 

Finally, we extrapolate to A; = cxd and obtain 

soo ~ 1.4999949. (7.4) 

This is very close to the expected result s = 3/2. The polynomial fit as well as the 
extrapolation is shown in Figure 6. 



0.002 0.004 0.006 0.008 0.01 




Figure 6. Polynomial fit to the exponents of (7.2) plotted against 1/k for 
a = f3 = 0.7. Extrapolation gives an exponent s^c = 1.5. 



7.2. Subleading corrections 

Having established that the leading behaviour of £i{L) at large L is as L~'^/^, we now 
turn to subleading corrections. Assuming that 

£i{L)^-eL-^/^-fL-'' , (L^«3), (7.5) 

we wish to determine the value of the exponent d. To this end, we define the sequence 
L 



(L + 2f/^£i{L + 2) - L^/^£i{L) 



(7.6) 
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If (7.5) is correct, we expect to be proportional to L^/'^^'^ at large L. In Figure 7 
we plot as a function of in a double-logarithmic plot. The result is well 
approximated by a straight line with slope 0.95, which suggests that the exponent of 
the subleading corrections is d = 5/2. 




Figure 7. Double logarithmic plot of as a function of 1/L for a = (3 ■ 
and 800 < L < 2600. 



0.7 



A more accurate estimate of the exponent can be obtained by extrapolating the 
finite-size data along the same lines as before. We group the data for into bins of 
11 points each and carry out a least-square fit for each bin to a functional form 



A, 



(7.7) 



The resulting sequence of exponents Sk is described very well by the polynomial least- 
square fit 



Sk « 0.998755255 - 4.329298118fc" 



.363104997fc" 



(7.8) 



Extrapolation gives Sqo ~ 0.998755255, which is very close to 1. This suggests that 
the subleading corrections to £i indeed scale like L^^/"^. The polynomial fit to the 
sequence of exponents Sk and the extrapolation to fc = oo are shown in Figure 8. 



0.004 0.006 0.008 0.01 



Figure 8. Polynomial fit to the exponents in (7.7) plotted against 1/k for 
a = = 0.7. 



7.3. Coefficient of the L term 

Having established that asymptotically the L-dependence of the eigenvalue is given 

by 

£i(L)--eL-3/2_/i-5/2 (L^oo), (7.9) 
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we now determine the coefHcient e. We again arrange the data for £i{L) into bins 
20fc < L < 20{k + 1) and within each bin perform least-square fits of £i{L) to the 
functional form 

« -efeL"^/^ - /fci"^/^ , 20fc < L < 20(fc + 1). (7.10) 

As is shown in Figure 9, the resulting sequence of coefficients for the largest available 
values of k (100 < k < 128) can be fitted very well to the polynomial 

Ck « 3.5780576 - 7.2704902fc-2_ (-j ll) 
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Figure 9. Polynomial fit to the coefficients ej. defined in (7.10) plotted against 
for a = /3 = 0.7. 

Extrapolation to fc oo then gives the following result for the eigenvalue of the first 
excited state of the TASEP in the MC phase 

£i{L) « -3.578 + ©(L^^/^). (7.12) 

As we have already mentioned, the behaviour (7.12) is in fact independent of a and 
P throughout the maximum current phase and the coefhcient is a universal number. 

The numerical value for the gap is smaller than that of the half-filled TASEP on 
a ring, where it is found that 5i,ring(i) 6.509 . ..L'^/^ [24,25,64]. 

7.3.1. Extrapolation Procedure In order to assess the accuracy of the numerical 
extrapolation procedure we have employed above, it is instructive to consider the 
analogous analysis on the coexistence line. Here the exact analytical value for the 
spectral gap is known. For a = /3 = 0.3 equation (6.18) gives 

£i(L) = -5.1815423 . . .L^^ + 0{L-^). (7.13) 

We have computed £i{L) on the coexistence line by a direct numerical solution of the 
Bethe ansatz equations for lattices of up to L = 1800 sites. Fitting £i{L) to the form 

£i(L) w -CfeL""''- , 20/c < L < 20(A:-|- 1), (7.14) 

and using the same analysis as above, but for the somewhat smaller values 61 < fc < 88, 
we find that the sequence of exponents Sk is well approximated by the polynomial 

sk ~ 2.0000967 - 0.3772271fc-^ + 0.7953103fc-2. (7.15) 

Extrapolation to A: = c» gives excellent agreement with the exact result s — 2. 
Similarly, a least-square fit of the sequence of coefficients Ck to a polynomial in 1/fc 
gives 

Cfc « 5.1814722 - 0.0002060fc~i + 0.0000782/;;-^ (7.16) 

The extrapolated value Coo agrees with the exact result (7.13) to five significant digits. 
This is quite satisfactory. 
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8. Other Gaps and Complex Eigenvalues 

Having established the large-L behaviour of the eigenvalue of M corresponding to the 
lowest excited state, we now turn to higher excited states. As M is not Hermitian, its 
eigenvalues arc in general complex. A complex eigenvalue in turn leads to interesting 
oscillatory behaviour in the time evolution. Aspects of such behaviour have been 
discussed for the KPZ equation [67, 68] 

8.1. Massive Phase I: a < ac, (3 < (3c, a ^ (3 

Here the next lowest excitation is obtained by choosing the integers Ij in the Bethe 
ansatz equations (5.6) as 

Ij = -L/2+j for j = l,...,L-2 /^.i = L/2. (8.1) 

This choice of integers Ij corresponds to a "hole" between the last two roots. Hence we 
will refer to this state as a "hole state" . As this choice is asymmetric with respect to 
the interchange j L—j, we expect the corresponding eigenvalue to be complex. This 
is indeed the case for small system sizes L. However, as the system size increases, the 
eigenvalue becomes real at a certain finite value of L. This comes about in the following 
way. For small L, the last root has a finite imaginary part and the distribution of 
the other roots is asymmetric with respect to the real axis, see Figure 10(a). As L 
increases the last root approaches the negative real axis until above a critical value of 
L its imaginary part vanishes. The other roots are then arranged in complex conjugate 
pairs, so that the corresponding eigenvalue becomes real. An example of this is shown 
in Figure 10(b). Further details regarding this phenomenon are given in Appendix C. 



• -0.5 2 4 

-2 

(a) • • • • 

Figure 10. Reciprocal root distributions corresponding to the hole state for 
a = 0.3 and /3 = 0.4; (a) L = 2n = 20 (b) L = 2n = 100. 

Assuming that the root distribution remains as in Figure 10(b) in the limit 
L — > oo, we can compute the corresponding eigenvalue in the following way. As 
Zl-1 is an isolated root, we first write the counting function in (5.6) as 

1 1 
iYL{z) = g{z) + -{g^iz) + K{zl-i, z)) + jY. ^(^'' (^-2) 

1=1 

We can now follow the same procedure as in Section 5.2 and turn this into an integro- 

diffcrcntial equation of the form (5.15). Importantly, the cndpoints of the integration 
contour are again complex conjugates of one another. Compared to the equation for 
the lowest excitation the driving term has an extra contribution of order 0{L^^) and 
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the boundary conditions that determine the endpoints ^ and ^* of the contour now 
read 



The root Q = zl-i is determined by the requirement that 

If we expand ^, ^* and Yl, but not (, in powers of L~^, we can use the results of 
Section 6. With 



(8.3) 
(8.4) 



OC 

YL{z) = J2L-^yn{z), 

n=0 



(8.5) 



we find the same solutions for yo, 2/2 and 2/3 as in (6.7), (6.10) and (6.11). However, 
due to the additional term proportional to in (8.2), the solution for yi differs from 
that in (6.9) and is given by 



iln 



z 1 



-iln 

Only now will we expand (, 



Zc 1 - 

z-l/C 



-iXi 



z + a z + b 



a Zc 



- 1/C. 



E 

n=l 



C = Co + > >-"<5n, 



+ M-iln(-r'a6(--2c)-'^0- (8-6) 



(8.7) 



and use the definitions of the points ^, ^* and ( above. Employing again the expansion 
(6.2) and Appendix B, we find that the intersect of the solution curve with the negative 
real axis in the limit i — > 00 is given, as before, by 

^ (8.8) 



ab 

Furthermore, in leading order, the isolated root C tends to Zc in the limit, i.e. Co = -^c- 
Putting everything together we finally find that the eigenvalue for this solution is given 

by 

S^iL) = -a-(3--^ + 0{L-% (8.9) 

with Zc given in (8.8). It is instructive to compare £2{L) to the gap of the first excited 
state £i{L). The latter exhibits a non-analytic change at /? = /3c. Interestingly, we 
find that to order 0(1.-2) 

lim £2iL) = lim £AL) 7^ lim £i(L) . (8.10) 

In conclusion we find that in the massive phase Mi also the second gap is real and 
given by (8.9). The second gap may be a complex conjugate pair for small values 
of the system size. However, this pair merges at a finite value of L producing two 
real eigenvalues, the lowest of which is given by (8.9). This observation is consistent 
with the results of Dudzinski and Schiitz [30], who computed the spectrum on the 
coexistence line a = /3 for small system sizes. 
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8.2. Coexistence line 

As we have seen above, the second gap of the TASEP is real for large L, even though 
it may be complex for small system sizes. On the coexistence line, the leading term 
of (8.9) vanishes and wo arc left with a second diffusive mode with an eigenvalue 
that vanishes as 0{L^^). A priori we expect a whole band of diffusive modes on the 
coexistence line, some of which should be given by the domain wall theory of [30], 



SniL) = — 



L2 



+ 0{L-% 



(8.11) 



with Zc given by (8.8). It is clearly of interest to know whether or not there are 
complex modes as well, which would result in oscillatory behaviour in the relaxational 
dynamics at large times. 

We observed numerically, that the first complex excitation that does not become 
real for a finite value of L in the massive phase Mi corresponds to the choice of integers 

Ij = -L/2 + j for i = l,...,L-3 

Il-2 = L/2 - 1, Il-1 = L/2, (8.12) 

i.e. there is one hole between the second and third last roots. The corresponding 
(reciprocal) root distribution again has an isolated root on the negative real axis, see 
Figure 11, which now lies inside the contour of integration. 




Figure 11. Reciprocal root distributions corresponding to the first complex 
eigenvalue £c{L) for a = (3 = 0.3 and L = 1860. 



Proceeding as in the preceding section, we write the counting function as 

1 1 
iYL{z) = g{z) + -(ffb(^) + K{zL-i, ^) + 7 E (^-13) 



1=1 



Turning the sum into an integral we arrive at the following integro-differential equation 

iYL{z) = g{z) + i {g^{z) + K{Ci,z) - if(C2, z)) 

+ ^ j^' K{w,z)Yi{w)dw 

1 / K{w,z)Y[{w) 1 f K{w,z)Y[{w) 

Here, Ci = ^L-i is the isolated root and ^2 corresponds to the position of the hole. The 
contribution due to the latter needs to be subtracted in order to cancel the contribution 
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arising from the integral. The values of Ci, C2 as well as the endpoints ^1 and ^2 of 
the curve are determined self-consistently by the "boundary conditions" 



YUi2) = -TT + -, 
Yl{Ci) = TT, 

Assuming an expansion of the form 

00 

71=0 



YL{C2)=7r- 



TT 

L' 

T' 



(8.15) 



(8.16) 



wc can utilize the results of Section 6. If we furthermore assume that C,^^ lies inside 
the contour of integration, we find that yo{z) is again given by (6.7) and obtain the 
following result for yi{z) 



yi{z) = - iln 



-iln 



^ 1 - 
Zcl- Z'^ 



Zc Zf. 



Z- Zc 



-iAi 



z-l/Ci z-Ci 



+ iln 



- 1/Cl Zc - Cl 

Ki -iln(-C2a6(-Zc)"'^0 



z + a z + b 
Zc + azc + b 



Z-I/C2 



I/C2, 



.17) 



It is important to note, that expression (8.17) has been derived under the assumption 
that z~^ lies outside the contour of integration+. Hence, we cannot use (8.17) to 
determine (i via (8.15). However, it turns out that (8.17) already contains sufficient 
information for determining the leading order of fc- 

Using the definitions of the points ^1, ^2 and (2 above, the expansion (6.2) and 
Appendix B, we find that the intersect of the solution curve with the negative real 
axis in the limit L ^ 00 is now given by 

Zc = -iab)-^/\ (8.18) 

The 0(1) term of the eigenvalue corresponding to the choice (8.12) is real 

1 -I- "^r 

Sc{L)=-a-P--^^+o{l), (8.19) 

I — Zc 

where Zc given in (8.18). The subleading corrections are complex. However, the 0(1) 
contribution (8.19) to the eigenvalue does not vanish on the coexistence line. Hence 
this excitation does not play an important role at large times. We have checked that 
for a = /? = 0.3 both (8.18) and (8.19) agree well with a direct numerical solution of 
the Bethe ansatz equations. Based on the above analysis, we conjecture that for all 
gapless excitations on the coexistence line with eigenvalues that scale with system size 
as 



£ 



-eL- 



0{L-^), (8.20) 

the coefficients e are real. In other words, the dominant large-time relaxation on the 
coexistence line does not have oscillatory modes. 

+ Otherwise the term ln(l — wz) in the kernel K(w,z) would produce an additional contribution in 
the integro-differential equation for yi{z). 
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8.3. Maximum Current phase 

Comparing again the numerical solutions of the Bethe ansatz equations to a direct 
numerical determination of the eigenvalues of the transition matrix for small system 
sizes L we find that the lowest excitations with non-vanishing imaginary part are 
characterized by the set of integers 

I^=-L/2+j for j = l,...,L-2 = L/2, (8.21) 

or 

Ij =-L/2+j for j = 2, . . . , L - 1 h = -L/2. (8.22) 

For the first choice the imaginary part of the eigenvalue is positive, whereas the second 
yields the complex conjugate eigenvalue. In contrast to the corresponding excitations 
in the massive phase Mi, we find that the gaps in the MC phase remain complex for 
all values of L. 

We denote the complex eigenvalue corresponding to (8.21) by £c{L). We have 
determined £c{L) from a numerical solution of the Bethe ansatz equations for lattices 
of up to L = 1200 sites and find its leading large L behaviour to be independent of the 
rates a and (3. In what follows we therefore fix q = /3 = 0.7, in the understanding that 
the results for the large-L asymptotics of the eigenvalue of M are universal throughout 
the maximum current phase. 

The distribution of reciprocal Bethe roots for L = 400 and a = /3 = 0.7 is shown 
in Figure 12. We observe that the root distribution is quite similar to that of the 
second lowest eigenvalue in Figure 3. The main difference is that now there is a slight 
gap between the last two roots. 




Figure 12. Reciprocal root distributions corresponding to the lowest complex 
eigenvalue £c(L) for a = /3 = 0.7 and L = 2n = 400. 

Based on the similarity of the root distribution to the one of the lowest excited 
state, we expect that both real and imaginary parts of £c{L) will be proportional 
to L~^/'^ for large systems. This is indeed the case, as is shown in Figure 13, where 
Re(fc(i)) and lTa{£c{L)) are plotted as ftmctions of L^^ on a double logarithmic scale 
for a = /3 = 0.7. Both curves are very close to straight lines with slope 3/2. 

In order to determine the asymptotic form of £c{L) more accurately, we repeat 
the analysis we employed for the lowest excited state in the maximum current phase 
in section 7. We group the data points for £c{L) in bins containing 11 points each and 
within each bin perform least square fits to 

Re{£^{L)) « Re(efe)L-«'= , 

Im(fc(i)) « Im(efc)L-*'' , 2Qk < L < 20{k + 1). (8.23) 
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Figure 13. Double logarithmic plots of the real and imaginary parts £c as 
functions of l/L. 

The resulting sequences of exponents Sk and tk in the range 30 < fc < 59 are well 
approximated by the polynomials 

Sk « 1.4997697 - 0.737888549fc-i + 3.369432984A:-2 , 

tk « 1.4991041 - 0.423377041fc"i + 6.133640280^-^. (8.24) 

The extrapolated values at A; = oo are very close to 3/2. The fits to the data as well 
as the extrapolation to A; = oo are show in Figure 14. 



1.502 




Figure 14. Fits to the exponents sj; and ti^ in (8.23) of the real and imaginary 
parts of £c{L) plotted against 1/fe for a = /3 = 0.7. 

In order to determine the coefficients of the 0{L^^^^) terms, we carry out least 
square fits of the binned data for £c{L) to 

fc(^) w „ ^ 20A: < L < 2G(fc + 1). (8.25) 

The resulting sequence of coefficients is well described by the polynomial fits (see 
Figure 15) 

Re(efe) « 8.4687424 - 0.1886671fc-i - 12.277915fc-2, 

Im(efe) w -1.6508588 + 0.1263865fc-i + 3.6333854A:-2. ^^'^^^ 

Extrapolation to the limit A: — > cx) then gives the following result for the energy of the 
first complex excited state of the TASEP in the MC phase 

£c{L) w -(8.47 - 1.65 i) L^^/^ ^ 0(^-5/2^^ 27) 

The result (8.27) can be compared to excited states with complex eigenvalues in the 
half-filled * TASEP on a ring. We are not aware of any explicit results in the literature 

* It is natural to compare the the half-filled case, as the average bulk density of the TASEP with 
open boundaries in the maximum current phase is 1/2. 
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Figure 15. Polynomial fits to the real and imaginary parts of the coefficient 
in (8.25) plotted against 1/fe for a = /3 = 0.7. 



on complex eigenvalues, but they can be easily determined by employing the method 
of [64]. We summarize some results in (Appendix D). The lowest excited state with 
complex eigenvalue found to be 

fc,ring(i) ~ -(17.1884 ... - 5.43662 . . . i) L'^^"^. (8.28) 



9. Summary and Conclusions 

In this work we have used Bethe's ansatz to diagonalize the transition matrix M for 
arbitrary values of the rates p, q, a, /3, 7 and 5 that characterize the most general 
PASEP with open boundaries. The resulting Bethe ansatz equations (3.1), (3.2) 
describe the complete excitation spectrum of M . 

We have carried out detailed analyses of the Bethe ansatz equations for the 
limiting cases of symmetric and totally asymmetric exclusion and determined the 
exact asymptotic behaviour of the spectral gap for large lattice lengths L. This gap 
determines the long time (t ^ L) dynamical behaviour of the TASEP. We emphasize 
that care has to be taken regarding time scales, and that our results below are 
not valid at intermediate times where the system behaves as for periodic boundary 
conditions [4]. 



9.1. Dynamical phase diagram 



In the case of totally asymmetric exclusion and 7 = 5 = 0, we found that there are 
three regions in parameter space where the spectral gap is finite and the stationary 
state is approached exponentially fast. In addition there is one region (maximum 
current phase) and a line (coexistence line) where the gap vanishes as i ^ cxd. The 
resulting dynamical phase diagram (j is shown in Figure 16. 
Recalling that 



1 



I- a 



1/3- 



1 



/3 



l-(3 



1/3' 



(9.1) 



the leading asymptotic values of the spectral gap in the various regions of the phase 
diagram of the TASEP are as follows: 



tt We use the term "phase" to indicate a region in parameter space characterised by a particular type 
of relaxational behaviour. 
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Figure 16. Dynamic piiase diagram determined by the first eigenvalue gap of the 
TASEP. Mj, Mjia a-iid ^llh are massive phases, CL denotes the critical coexistence 
line and MC the critical maximal current phase. 



Massive Phase Mi: a < ac, (3 < (3c, a ^ (3 

£,{L) = -a-f3+ \ == + 0{L-\ (9.2) 

1 + ^(1 - a)(l - l3)/a[3 

The spectral gap does not vanish as L — > oo, which impUes a finite correlation length 
and an exponentially fast approach to stationarity. 

High-Density Phase Mub'. (3 < 1/2, ac < ot 

£i{L)^-ac- P+ ^ — + 0(L-2). (9.3) 

I + [{I -13) IP]'" 

The spectral gap is finite and independent of [3. The relaxational behaviour is again 
exponentially fast. 

The gap in the low-density phase Mnai a<l/2, /3c</?is obtained by the 
exchange a /3. We note that the subdivision of the massive high and low density 
phases into Mi and Mna,b is different from the one suggested on the basis of stationary 
state properties in [16]. 

Coexistence Line (CL): (3 — a < 1/2. 

f,(L).-^^Mi-i^L- + 0(L-). (9.4) 

1 — 2q; 

The gap vanishes like for large systems, which corresponds to a dynamic exponent 
z = 2 and indicates relaxational behaviour of a diffusive type. 

Maximal Current Phase (MC): a,(3 > 1/2. 

£i{L) w ^ 0(^-5/2^^ (9 5^ 
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Here the gap is independent of the rates a and /3. The dynamic exponent z = 3/2 is 
indicative of KPZ like behaviour [5]. Wc find that the magnitude of £i{L) is smaller 
than half that of the lowest excited state for the TASEP with periodic boundary 
conditions 

^i,ring(i) ~ -6.50919 . . . (9 g) 

The two gaps do not seem to be related in any obvious way. 

It is known [17] that by varying the bulk hopping rates one can induce a 
crossover between a diffusive Edwards- Wilkinson (EW) scaling regime [18] with 
dynamic exponent z = 2 and a KPZ regime [5] with z = 3/2. Here we have shown 
using exact methods that a crossover between phases with z = 2 and z = 3/2 occurs 
in the case where the bulk transition rates are kept constant, but the boundary 
injection/extraction rates are varied. 



9.2. Domain wall theory 

As shown in [29,30] the diffusive relaxation {z = 2) is of a different nature than in 
the EW regime and is in fact due to the unbiased random walk behaviour of a shock 
(domain wall between a low and high density region) with right and left hopping rates 
given by 

P+ - P 

Here, p~ — a and /9+ = 1-/3 arc the stationary bulk densities in the low- 
and high-density phases respectively, and the corresponding currents are given by 
p^{l — p^). Our results (9.2) and (9.4) for the massive phase Mi and the 
coexistence line agree with the relaxation time calculated in the framework of a domain 
wall theory (DWT) model in [30]. This is in contrast to the massive phases Mu, where 
the exact result (9.3) differs from the DWT prediction of [30]. More precisely, the 
DWT predicts that the relaxational mechanism in the stationary high-density phase 
(/? < 1/2, a > (3) for a < 1/2 is due to the random walk of a "(G|l)" domain wall 
between a low and a high density segment. On the other hand, in the high-density 
phase for a > 1/2, the DWT predicts a relaxational mechanism based on so-called 
"(m|l)" domain walls between a maximum current and a low-density region. Our 
results for the gap exhibit a change of behaviour at a = ac rather than at a = 1/2. 
We therefore propose the following modified DWT. 

We assume that in the high density phase for large a we can still use the DWT 
rates as given in (9.7), but with an effective density p~g. In contrast to [30] we do 
not take the effective density equal to that of the maximum current phase, = 1/2, 
but instead determine it below from using a monotonicity argument. Consider first 
Figure 17, in which we plot the first gap as a function of p~ = a for constant /3 = 0.3. 

• a < /3 

We are in the low density phase. The gap in the infinite volume is finite in this 
region, but upon increasing a we are driven towards the coexistence line where 
the gap vanishes and the relaxation is diffusive. Hence, d£/da > in this region. 

• a > /3 

We are now in the high density region and we expect the gap to be finite again, 
decreasing with increasing a. Hence, in this region we should have d£/da < 0. 
Initially, this is indeed the observed behaviour of the graph of (9.2), but we 
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see that at a — the graph of (9.2) has local minimum, where the expected 
behaviour breaks down. 




Figure 17. The gap as a function of a for /3 = 0.3. The curve corresponds to the 
function in (9.2) and the horizontal line to (9.3). The gap of the lowest excitation 
(red online) is a combination of the curve (0 < a < etc) and the line (qc < a < 1). 



Hence, for small values of p and large values of p+ the domain wall theory of [30] 
correctly predicts the gap to be 



£i{p-,p+) = -D+ - D- +2VD+D- , (9.8) 

which is equal to (9.2). In the high-density region, there is a critical value 
beyond which the gap is given by £i(/3^jj, p^), where p~g is determined by 

d£i{p-,p^ 



dp 



= 0. (9.9) 

The considerations above are of a sufficiently general nature to remain valid for the 
wider class of driven diffusive systems that can be described by effective domain wall 
theories. We emphasize that the DWT as discussed above pertains to the limit L — > oo 
only, and that a more careful analysis is necessary to correctly produce the difference in 
finite size behaviour for the phases Mi and Mu, see (6.15) and (6.23). In particular, 
the boundary rates generally need to be chosen differently from the bulk rates. A 
more detailed comparison of the Bethe ansatz solution and DWT is possible for the 
case of the PASEP with Q ^ 0, where DWT becomes exact for certain values of the 
rates [65] . This is beyond the scope of the present work and will be discussed in a 
separate publication [66]. 



9.3. Higher excitations 

Massive Phase Mi and Coexistence Line (CL): 
The second gap in this phase is given by 

2zr 47r2 



£2{L) 



-a- 13- 



1 



L2 



0{L- 



(9.10) 
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where Zc = —1/Vab. This result is in agreement with the second gap predicted by 

the DWT of [30]. In particular, on the coexistence line it corresponds to n = 2 of the 
DWT band of diffusive modes on the coexistence line, 

£n{L) = -^^-^ + 0{L-^). (9.11) 

Zc — Zc ^ 

The first gap given in (9.4) corresponds to n = 1. 

We studied the possibility of an oscillating diffusive mode on the coexistence line 
and found that the first excitation with nonzero imaginary part has a finite gap for 
L — > oo 

£,{L) = -a-P-^-^+oil), (9.12) 

1 Zc 

where Zc = —{ab)~^^^. Although the subleading terms have nonzero imaginary parts, 
the leading term is real and does not vanish on the coexistence line. Hence this 
excitation does not play an important role at large times. We conjecture that the 
leading 0{L~'^) terms of the eigenvalues of all diffusive modes on the coexistence line 
are real. 



Maximum Current Phase MC 

We determined the energy of the first "complex" excited state in the MC phase by 
means of a direct numerical solution of the Bethe ansatz equations. We find it to be 
independent of a and /3 and given by 

£c{L) w -(8.47 - 1.65 i) L'^'"^ + 0{L-^'^). (9.13) 

Complex eigenvalues give rise to interesting oscillations in correlation functions. Such 
behaviour has been observed for the KPZ equation [67,68], which is related to the 
TASEP with periodic boundary conditions. In Appendix D we have computed the 
low lying complex excitations of the TASEP on the ring using the method of [64] . The 
lowest such excitation is given by 

^c,ring(i) ~ -(17.1884 ... - 5.43662 . . . i)L-^/'^. (9.14) 

The result (9.14) confirms the prediction from the KPZ equation that both real and 

imaginary part scale with L"'^/^ [68]. According to (9.13) this scaling still holds when 
boundaries are present. However, as was the case with the first gap, although (9.13) 
does not depend on the boundary rates a and /3, its value does not seem to be related 
to (9.14) in a simple way. 

A number of interesting open problems remain. In a forthcoming publication [66] 

we determine the spectral gaps for the case of partially asymmetric diffiision and 
Q,/3, 7,(5 from the Bethe ansatz equations (3.2). As we have discussed, the 
constraint (1.8) precludes the determination of current fluctuations from the Bethe 
ansatz. It therefore would be highly desirable to obtain Bethe ansatz equations for the 
Heisenberg XXZ chain with the most general open boundary conditions. Furthermore, 
the existing Bethe ansatz solution provides information only about the spectrum, 
but not the eigenstates of the transition matrix. In order to study correlations 
functions [39] the knowledge of matrix elements of the spin operators between left 
and right eigenstates is required. In light of this it would be very interesting to 
construct the eigenstates from the Bethe ansatz. 
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Appendix A. Analysis of the Abel-Plana Formula 



In this appendix we sketch how to extract the finite-size correction terms from the 
integral expression 



iYLiz) = g{z) + ^g^{z) + ^ J^^ K{w, z)Yl{w) dw 



+ 



1 r R 



K{w,z)Y[{w) 



K{w,z)Y[{w) 

(ULYi(«,) 



..^ , , , ■ , > dw. (A.l) 

q-iLYl{w) 27r 7^2 o'^^i("') - 1 

The main contributions to the correction terms in (A.l) comes from the vicinities of the 
endpoints Along the contour Ci the imaginary part of the counting function is 

positive, lin{YL{w)) > 0, whereas along the contour C2 it is negative, lin{YL{w)) < 0. 
As a result the integrands decay exponentially with respect to the distance from the 
endpoints. In the vicinity of we therefore expand 

Yl{w) = Yl{0 + YliOiw - + • • • • (A.2) 
Assuming that Y^(^) is 0{1) (an assumption that will be checked self-consistently), 
we find that the leading contribution for large L is given by 

J_ f K{w,z)Y[{w) Y[{0 K{w,z) 

(A.3) 



dw 



Here we have used the boundary conditions (5.14). Carrying out the analogous analysis 
for the integral along the contour C2, we arrive at the following expression for the 
leading contribution of the last two terms in (A.l) 
J_ f K{w,z)Yl{w) 



A 



dw + 



K{w,z)Yaw) 



1 



(^ 



K £ + 



LYm' 



dw 



K 



LY[iO' 



dx 

(A.4) 



2nL Jo 1 

If the endpoints ^,^* are such that we can Taylor-expand the kernels appearing in 
(A.4), we can simplify the expression further with the result 

K'{^,z) r X 



f 

Jo 



^L^YLiO JO 

TT K'itz) 



1 + e^ 

- (e- 



dx - (^ ^ r) 



(A.5) 



12^2 YliO 

This is the leading Euler-MacLaurin correction term that occurs in the low and high 
density phases, see (6.1). The key in the above derivation was the ability to expand 

ix \ ^/^ ix 



K 



In 



LYliOiz- 







1 X 



LYliOiz-'-O-ix 



(A.6) 
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in a power series in x. This is unproblematic as long as Z/Y^(^)(z~^ — ^) is large, 
which turns out to be the case in the low and high density phases as well as on the 
coexistence line. 

In the maximum current phase the above analysis no longer holds. As mentioned 
in the main text, ^ ^ — 1 in the MC phase, and this has two consequences. Firstly, 
by virtue of the endpoints ^,^* being close to stationary point, the derivative Yl{£^) is 
now found to be of order L~^/^ whereas Y^{(,) is of order one. As a result (c.f. [25]) we 
have to retain subleading terms in the Taylor-expansions of Yl{w) and Y^{w) around 
w = ^, e.g. 

yl{w) ^ yl{o + yaoi^ -0 + ^Y^iOi^ - o'- (a.7) 

Secondly, and more seriously, as ^ — = 0{L^^^) the Taylor expansion of the 
integration kernel K{^, z) breaks down near z = Taking (A.7) into account and 
following through the same steps as before, we find that the leading contributions of 
the last two terms in (A.l) in the maximum current phase are 

1 / X + p 



A=— kU 

ttL J a 1 + exp(a;2 + 2xp) V 

X — ip 



+ exp(a;2 — 2\px) 
Here we have introduced 



K(^^ + iex,z) - (^^f).(A.8) 



0{L-"^), p=-^sY[{0L = O{l). (A.9) 



In order to make progress we would like to expand around x = 0. However, as — ^ 
becomes of order 0{L^^^'^) for z near ^, all the terms x"i\r("^(^, z) in the expansion 
are of the same order L'^ for [n > 1), and we would need to self-consistently determine 
the full series. The n = terms in the Taylor expansion, which are proportional to 
ln(L)L^^, cancel, and we find that A — 0{L^^). 

Appendix B. Expansion coefficients 

In the following list we abbreviate y'nizc) by y'^, 

Ki=-y'oSu (B.l) 
^i = yo-, (B.2) 

«2 = - y'oS2 - y[di - Ivoi^l - Vl), (B.3) 

^2 = - {y'om + y'lVi + y'oSim) , (b.4) 

TT 

M2 = 2/^^, (B.5) 

TT 

K3 = - y'oh - v'lh - y'25i - 2/o (<5i^2 - mm) - \yii^i - vl) 

- Jyo%(^i-37?D, (B.6) 

As = - {y'om + y'lm + y'2'ni + y'oi^im + ^2??i) + ^imy'l 

TT 
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1 



If 1 

= - y'oihm + <52??i) + y[5ir]i + -2/or/i(3^f - rii) 



+ —772^1 



^-3 = 



3 y'^ 

1 (.^..(3.?-.?) -2.^1). 



(B.7) 

(B.8) 
(B.9) 



Appendix C. Complex/Real excited states and their Bethe root 
distributions 

As mentioned in the main text, excited states can change from complex- to rcal-vahied 
with increasing lattice length. In order to illustrate this point we consider the TASEP 
with a = (3 = 0.25. In Fig. CI we plot the real parts of the eigenvalues of the transition 
matrix for lattice lengths L = 4, 5, . . . , 9, computed by direct diagonalisation. The 
lowest excitation always corresponds to a real eigenvalue. For X = 4, 5, 6 the excited 
states with second largest real part are a pair of complex conjugated eigenvalues. 
For larger L however, this pair becomes a pair of real eigenvalues. In terms of the 
Bethe roots this comes about in the following way. The pair of complex conjugated 
eigenvalues corresponds to distributions of integers 

Ij = -L/2 + j for j = l,...,L-2 = L/2, (C.l) 

and 

for 



h - ~L/2, 



(C.2) 



I,=^L/2 + j for J = 2,. 

respectively. The numerical values of the corresponding Bethe roots for even L are 
listed in Table CI. When the lattice length exceeds L = 6 the nature of the excited 
states changes. Now there are two real eigenvalues corresponding to the distribution 
(C.2) of integers, see Table CI for the L = 8 site system. Analogous violations 
of the one-to-one correspondence between sets of integers and roots of the Bethe 
ansatz equations have been previously observed in both "ferromagnetic" [69] and 
"antiferromagnetic" [70] situations. 




Figure CI. Real parts of the eigenvalues of low-lying excited states for small 
system sizes L = 4,5,. ..,9. The lines arc a guide to the eye only. Dashed lines 
indicate complex conjugate pairs of eigenvalues, solid lines correspond to real 
eigenavlues. 
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Zz 


= 


0.193853-0.0505861 






1 


Zi 


= 


0.190087+0.0724531 






2 


Zb 


= 


0.130508 + 0.2046321 


8 


-0.156354 


4 


Zl 


= 


-0.549451 






2 


Z2 


= 


0.062316 + 0.2291561 






1 


Zz 


= 


0.141563 + 0.1346571 









za 




0.173073 + 0.0442591 






-1 


Zb 


= 


0.173073-0.0442591 






-2 


zq 




0.141563- 0.1346571 






-3 


Z7 




0.062316-0.2291561 


8 


-0.218262 


4 


Zl 




-1.20057 






2 


Z2 




0.0876066 + 0.2527701 






1 


Zz 




0.163307+0.1447761 









Zi 




0.193253 + 0.04716571 






-1 


Zb 




0.193253 - 0.04716571 






-2 


Zq 




0.163307-0.1447761 






-3 


Z7 




0.0876066 - 0.2527701 



Table CI. Low-lying excited states for a = f} = 0.25 and L = 4, 6, 8. 



Appendix D. Excited States of the TASEP on a Ring 

In this appendix we collect some results on the excitation spectrum of the TASEP 
on a ring, I.e. with periodic boundary conditions. Gollnelll and Malllck have recently 
developed a simple method for calculating the spectral gap in the TASEP on a ring [64] . 
Their method applies to higher excited states as well. In what follows we consider the 
case of half-filling only, i.e. L/2 particles on a ring with L sites. The Bethe ansatz 
equations of the half-filled TASEP are [24, 25] 

(l-Z|)V^ = -2-nt^' ^■ = l> (D-l) 

"-'^ Z 1 
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As the right hand side of (D.l) is independent of the index j, we may write [24,25,64] 

(1 - ^ _g7r« ^ _^ < jjj^(y) ^ ^ p 3) 

The L roots of this equation are 

Zm = -ZL/2+m = \/l " 2/m , m=l,...,^, (D.4) 

where 

The j/m's he on a circle with radius e^'^"/^ such that 

< arg(yi) < arg(y2) < • • • < arg(yL/2) < 27r . (D.6) 

In order to construct a particular solution of the Bethe ansatz equations, one chooses 
a sequence of roots Vdi), ■ ■ ■ tVc{l/2) where 1 < c(l) < ... < c(n) < L and then 
determines the parameter u self-consistently from 

e-» = 2^ n . (D.7) 



Appendix D.l. Ground State 

The ground state of the TASEP on a ring is obtained by choosing c(j) = j, 
j = 1,. . . ,L/2. The corresponding distribution of roots ym is shown in Figure Dl. 




Figure Dl. Distribution of roots ym corresponding to the ground state of the 
TASEP on a ring. 



Appendix D.2. Particle-Hole Excitations 

Some simple excited states can be constructed by choosing sequences of the type 

{j if J < m 

3 + 1 if m < j < f (D.8) 
L-k iij = ^. 

Such excitations correspond to having a hole at ym and an extra particle at yL-k as 
compared to the ground state. The roots distribution is shown in Figure D2(a). In 
what follows we are interested in the limit i — > oo, while keeping k and m fixed. 

Golinelli and Mallick [64] have shown how to simplify and solve the self- 
consistency condition (D.7) in the large-L limit for the first excited state, which 
corresponds to the choice (m, fc) = (1,0) in (D.8). Following through exactly the 
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Figure D2. Particle-hole excitations with energies (a) and (b) 

EL/2-m.,L/2+k- 



same steps for general {m,k), we obtain the following self- consistency condition for 
the parameter u in the limit L — > oo 



Li3/2(-R''") +27r v^-u - (2m - l)i + y^-u + {2k + l)i = 0. (D.9) 

Once u is known from a numerical solution of (D.9) the leading behaviour of the 
eigenvalue of the transition matrix can be determined from 



Em,L-k = 0'm,L-k{u) 



3/2 



+ 



(D.IO) 



^m.L — k 



{u) 



Jill 



{-u - (2m - l)i)^^^ + {-u + (2fc + l)i) 



3/2 



1 



Li5/2(-e"" 



)• 



The momentum of the particle-hole excitations is given by 



-{m + k) — . 
Ij 



(D.ll) 



(D.12) 



We have solved (D.9) for the first few excited states of this type and list the 
corresponding values of u and the coefficients am,L-k that characterize the asymptotic 
behaviour eigenvalues of the transition matrix in Table Dl. 

The choices (to, k) = (2, 0) and (m, k) = (1, 1) give the eigenvalues with maximum 
real part and non-zero imaginary parts. 

A second sequence of simple particle-hole excitations is obtained by the choice of 
sequence 



c{3) 



3 

J + 1 



if.?<f 

if # - TO < J < f 



2 

m 

f+fc ifi = f. 



(D.13) 



This corresponds to having a hole at yL/2-m an extra particle at yL/2+k- It is 
easy to see that energy and momentum of such excitations are given by 

EL/2-m,L/2+k = ■E'fc,L-m> (D-14) 

PL/2-m,L/2+k = {m + k) — . (D.15) 
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(m, k) 


u 


^m,L — k 


(1,0) 


1.11907 


-6.50919 


(2,0) 


1.55661 - 0.18053i 


-17.1884- 5.43662i 


(1,1) 


1.55661 + 0.180531 


-17.1884 + 5.436621 


(2,1) 


1.90545 


-28.9435 


(1,2) 


1.84895 + 0.313567i 


-30.6202 + 14.05231 


(3,0) 


1.84895 - 0.313567i 


-30.6202 - 14.05231 


(2,2) 


2.15672 + 0.1414191 


-43.0859 + 8.573381 


(3,1) 


2.15672 - 0.1414191 


-43.0859- 8.573381 


(1,3) 


2.07629 + 0.4221661 


-46.2912 + 25.16771 


(4,0) 


2.07629 - 0.422166i 


-46.2912- 25.16771 



Table Dl. Eigenvalues of some low-lying excited states for the half-filled TASEP 
on a ring. 



Appendix D.3. Multiple Particle-Hole Excitations 



Multiple partlcle-holc excitations can be constructed along the same lines. We note 
that while the momentum is additive, i.e. simple the sum of the momenta of the 
constituent particle-hole excitations, this Is not the case for the energy. Hence the 
Bethc-ansatz particle and holes still Interact with one another. Let us first consider a 
two-particle two-hole excitation characterized by the sequence 



c(j) = { 



3 

J + 1 
i + 2 
f +fc' 
L-k 



If j < m 
if TO < j < 



if ^ 



m < i < f - 1 



This corresponds to having two holes at positions and 
particles at yL/2+k' and yL-k, see Figure Appendix D.3(a). 



yL/2- 



(D.16) 



and two extra 



L/2+k 
U2+W 




(b) (c) 
Figure D3. Two-particle two-hole excitations. 



Following once more the procedure of [64] we arrive at the following equation 
determining the parameter u in the L oo limit 



Li 



3/2 ( 



2n 



v/-it- (2to- 1)1+ + (2A: + 1)1 

= 0. 



+ y/-U + (2to' + 1)1 + y/-U - {2k' - 1)1 



(D.17) 
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(to, k, to', k') 


u 




(1,0,0,1) 


1.65874 


-16.0176 


(1,1,0,1) 


1.99633 + 0.124903i 


-27.8453 + 5.011651 


(2,0,0,1) 


1.99633- 0.124903i 


-27.8453-5.011651 


(1,0,1,1) 


1.99633 + 0.1249031 


-27.8453 + 5.011651 


(1,0,0,2) 


1.99633- 0. 1249031 


-27.8453 - 5.011651 


(2,0,0,2) 


2.31089 - 0.2055441 


-40.3126 - 9.523281 


(1,1,1,1) 


2.31089 + 0.2055441 


-40.3126 + 9.523281 


(2,0,1,1) 


2.29142 


-40.4698 


(1, 1,0. 2) 


2.29142 


-40.4698 


(1.0, 1. 2) 


2.29142 


-40.4()9(S 


(2,1,0,1) 


2.29142 


-40.4698 


(1,2,0,1) 


2.23876 + 0.2398581 


-42.0218+ 13.32191 


(1,0,2,1) 


2.23876 + 0.2398581 


-42.0218+ 13.32191 


(2,1,1,1) 


2.57296 + 0.08815721 


-53.6251 + 4.56691 


(1,1,1,2) 


2.57296 + 0.08815721 


-53.6251 + 4.56691 


(2,0,1,2) 


2.57296- 0.08815721 


-53.6251 - 4.56691 


(2,1,0,2) 


2.57296 - 0.08815721 


-53.6251 - 4.56691 



Table D2. Eigenvalues of some low-lying two-particle two-hole excited states for 
the half-filled TASEP on a ring. 



The transition matrix eigenvalue of the two-partlcle two-hole excitation for very large 
L Is given by 

1 3/2 



Em,L-k,L/2-m',L/2+k' — ('"m,L-k,L/2-m',L/2+k'{u) 



+ ...(D.18) 



^3/2 , , 
am,L-k,L/2-m',L/2+k'{u) = -^|(~"~ (2to - l)l) 

+ {-U + {2k + 1)1)'/' + (-« + (2to' + 1)1)'/' 
+ (_„ _ (2fc' - l)i)3/2} - _i_Ll5/,(-e-). 

The momentum Is equal to 



Pm,L-k,L/2-m',L/2+k' = {m' + k' - m - k) 



277 



(D.19) 



(D.20) 



We have solved (D.17) and (D.19) for some low-lying excited states and list the results 
In Table D2. 

The two-partlcle two-hole excitations shown In Figures Appendix D.3(b) and (c) 
can be c;onstmcted along the same lines. The transition matrix eigenvalues In the 
large-L limit fulfil 



-^L/2-m,L/2+fe,L/2-m',L/2+/c' = E, 



m,L-k,L/2-k' ,L/2+m' , 
'k,L-m,L/2-m',L/2+k' ■ (D.21) 
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